Accepted for publication in the Spitzer Edition of the Astrophyical Journal Supplement Series 

Preprint typeset using I^T^X style cmulatcapj v. 2/19/04 



SPITZER-IRS SPECTROSCOPY OF THE PROTOTYPE WOLF-RAYET STAR EZ CMA (HD 50896) 

Patrick W. Morris 1 

Spitzer Science Center, IPAC, California Institute of Technology, M/S 220-6, 1200 E. California Blvd., Pasadena CA 91125 

pmorris@ipac . caltech . edu 



o 
o 

(N 

c 



> 

o 

O 

o 

9*'- 
6 



X 



Paul A. Crowther 

Department of Physics & Astronomy, University of Sheffield, Hicks Building, Hounsficld Rd, Sheffield S3 7RH, United Kingdom 

AND 

Jim R. Houck 

Astronomy Department, Cornell University, 106 Space Sciences Bldg., Ithaca, NY 14853 
Accepted for publication in the Spitzer Edition of the Astrophyical Journal Supplement Series 

ABSTRACT 

We present mid-infrared Spitzer-IRS 2 spectroscopy of the prototype WN star EZ CMa (HD 50896, 
WN4b). Numerous stellar wind lines of Hell are revealed, plus fine-structure lines of [Nem] 15.5/im 
and [O iv] 25.9/im. We carry out a spectroscopic analysis of HD 50896 allowing for line blanketing and 
clumping, which is compared to the mid-IR observations. We make use of these stellar properties to 
accurately derive Ne/He=1.2-1.8xl0 -4 and O/Hc=4.2-4.8xl0~ 5 by number, for the first time in an 
early WN star. In addition, we obtain N/C^40 and N/O~50 by number, values in perfect agreement 
with current predictions for rotating massive stars at the end of interior hydrogen burning. 
Subject headings: techniques: spectroscopic — stars: individual (HD 50896) — stars: atmospheres — 
infrared: stars 



1. INTRODUCTION 

Wolf-Rayet (WR) stars are descended from the most 
massive O stars (Mi n ; t > 25M ), and shed core- 
processed material in thick, line-driven winds at such 
prodigious rates (typically 10~ 5 — 1O~ 4 M yr _1 ) as to af- 
fect the evolution of the core and thus the lifetime of the 
star, prior to ending in a supernova explosion (Meynet 
& Maeder 2003). The prototype Wolf-Rayet star EZ 
CMa (HD 50896 = WR 6 in the catalogue of van der 
Hucht 2001) has been studied extensively at nearly all 
wavelengths accessible with ground-based facilities and 
space observatories, to understand the extreme stellar 
wind properties of these rare massive stars, and how they 
interact with their local environment. 

At infrared wavelengths, the continuum opacity is prin- 
cipally determined by free-free scattering of electrons 
from helium ions, and the emergent radiation originates 
in the outer layers of the wind where the asymptotic ve- 
locity is reached (e.g., Hillier et al. 1983). These layers 
are also where critical electron densities and tempera- 
tures are reached to form fine structure lines of Ne, O, S, 
and other elements, that are instrumental in testing evo- 
lutionary model predictions of the surface abundances 
during different core burning stages. 

Elemental abundances obtained from the fine structure 
lines depend on the stellar parameters, particularly the 
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mass-loss rat, which is sensitive to the degree of clumping 
in the wind. The role played by mid-infrared spectra, in 
conjunction with UV, optical, and near-IR spectroscopy, 
has been demonstrated for a few WN and WC stars ob- 
served with the SWS (Morris et al. 2000; Dessart et al. 
2000). Only the brightest WR stars could be observed 
with the SWS, for sensitivity reasons, and observations 
of only five WR stars produced spectra of suitable quality 
for modelling beyond 4.5 /im. 

In this Letter we present Spitzer-IRS spectra of the pro- 
totype early WN star EZ CMa. We carry out an analysis 
of the ultraviolet to mid-IR spectroscopy, revealing stel- 
lar and wind parameters, with which Ne/He and O/He 
abundance ratios are determined. Comparisons with the- 
oretical expectations are presented. 

2. OBSERVATIONS 

EZ CMa was observed with the Spitzer InfraRed Spec- 
trometer (described by Houck ct al. 2004) on 11 Novem- 
ber 2003, during the Science Verification phase, using all 
four of the IRS modules. The data were acquired with 
the Spectral Mapping AOT, with the star observed at 3 
to 5 discrete positions along the lengths (cross-dispersed 
axes) of the slits. The 10 — 37.5 fim data presented 
in this paper are from the Short and Long High (SH 
and LH) resolution modules (R = A/AA ~ 700), and 
5.3 — 10 /im data are from the Short Low (SL) module 
(R ~ 70—120), using exposure times that should produce 
continuum signal-to-noise (S/N) ratios of around 100 in 
the high resolution spectra, and 200-250 in the low reso- 
lution spectra. The basic calibrations were automatically 
carried out in the SSC pipeline, version S9.5, and then 
spectra were extracted and coadded interactively, using 
the offline pipeline. 

Our analysis also makes use of archival Infrared 
Space Observatory (ISO) Short Wavelength Spectrom- 
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to (15-14) at 34.6/im. In addition to these, we observe 
fine-structure lines of [Oiv] 25.9/zm and [Nein] 15.5/im, 
the latter blended with Hell 15.47/xm (15-13). There 
is no evidence for [Siv] 10.5^m in HD 50896. The ob- 
served feature is well reproduced by a weak blend of He n 
10.46/im (21-15) and 10.50/im (24-16) in our synthetic 
spectrum. 
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Fig. 1.— Mid-IR spectroscopy of EZ CMa (HD 50896) from 
ISO/SWS and Spitzer/IKS with the principal lines identified. 



eter (SWS) spectra, optical spectrophotometry from 
Torres-Dodgen & Massey (1988), high resolution In- 
ternational Ultraviolet Explorer (IUE) ultraviolet spec- 
trophotometry from Howarth & Phillips (1986), plus 
Hopkins Ultraviolet Telescope (HUT) far-ultraviolet 
spectrophotometry (Schulte-Ladbeck et al. 1995). The 
unpublished SWS 2.4-4.2 fim spectroscopy of EZ CMa 
was obtained from the public archive, processed to the 
Auto Analysis Result. Further reduction to a single spec- 
trum was carried out interactively, using SWS Interactive 
Analysis routines. Continuum S/N ratios are around 50 
in these data. Low resolution near-IR spectroscopy of 
HD 50896 has previously been discussed by Hillicr et al. 
(1983). We additionally use high dispersion spectroscopy 
of HD 50896 obtained with the ESO 3.6m cassegrain 
echelle spectrograph (CASPEC), previously presented by 
Schmutz (1997). We adopted a distance of 1.8 kpc to 
HD 50896 (Howarth & Schmutz 1995), interstellar red- 
dening of E(B-V)=0.1 mag (the mean from Schmutz & 
Vacca 1991, van der Hucht 2001) and an atomic hydro- 
gen column density of \og(NHi) =20.7 cm -2 (Howarth 
& Phillips 1986). 

The Spitzer IRS and ISO SWS spectroscopy of EZ CMa 
are presented in Fig. 1 together with identifications of the 
principal wind features. The majority of mid-IR spec- 
tral features are due to wind lines of Hen, primarily the 
leading member in each series from Hen (7-6) at 3.09/im 



3. SPECTROSCOPIC ANALYSIS OF HD 50896 
3.1. Technique 

To date, the majority of spectroscopic analyses of WN 
stars have been carried out using non-LTE models that 
did not account for metal line blanketing or clumping 
(e.g. Crowther et al. 1995, Hamann et al. 1995). In 
only a few cases have clumped, metal line blanketed mod- 
els been applied to individual stars (e.g. Schmutz 1997, 
Morris et al. 2000, Dessart et al. 2000, Herald et al. 
2001). Such effects need to be taken into account in 
order to determine the fundamental parameters of WR 
stars (Crowther 2003). 

For the present study of HD 50896, we employ CM- 
FGEN (Hillier & Miller 1998), which solves the trans- 
fer equation in the co-moving frame subject to statisti- 
cal and radiative equilibrium, assuming an expanding, 
spherically-symmetric, homogeneous and static atmo- 
sphere, allowing for line blanketing and clumping. The 
stellar radius (i?*) is defined as the inner boundary of 
the model atmosphere and is located at Rosseland optical 
depth of ^20 with the stellar temperature (T*) defined by 
the usual Stcfan-Boltzmann relation. Consequently, the 
stellar radius is much smaller than the radius of t=1, 
such that the temperature depends directly on the as- 
sumed velocity law. 

Our approach follows previous studies (e.g. Crowther 
et al. 1995), such that diagnostic optical lines of Hci 
(A5876), Hen (A4686) plus the local continuum level al- 
low a determination of the stellar temperature, mass-loss 
rate and luminosity. We were unable to employ solely 
mid-IR diagnostics since the IRS spectrum of HD 50896 
is dominated by lines of He II, with He I weak or blended. 

Our final model atom contains H, He, C, N, O, Ne, 
Mg, Si, P, S, CI, Ar, Ca, Cr, Mn, Fe and Ni. In to- 
tal, 1186 super levels, 4462 full levels and 45,485 non- 
LTE transitions are simultaneously considered. We as- 
sume hydrogen is absent, such that helium makes up 
98% of the atmosphere. CNO abundances are varied 
until an optimal fit is achieved. In contrast with nitro- 
gen and carbon, the abundance of oxygen in WN stars 
has proven to be rather challenging (Hillier 1988). The 
reason is that oxygen has a much more complicated ion- 
ization stratification in their inner winds, plus diagnos- 
tic lines are located in rather observationally inaccessible 
regions, typically 2800-3400A. Consequently we set oxy- 
gen (and neon) abundances from the fine-structure anal- 
ysis. Other elements arc fixed at Solar values (Grevesse 
& Sauval 1998; Asplund et al. 2004). 

We adopt a standard form of the velocity law, v(r) = 
Woo(l — Ri,/r)P, where (3=1. In contrast, Schmutz (1997) 
solved the hydrodynamical equation for outer wind and 
obtained [3 ~ 3 with reference to the core radius. A 
terminal wind velocity, v^, of 1860 kms -1 is obtained 
from the mid-IR fine structure line of [O iv] 25.9/im. For 
comparison, Prinja ct al. (1990) obtained 1720 kms -1 
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from UV observations of the P Cygni Civ 1550A line, 
whilst Schmutz (1997) derived 2060 kms -1 from a fit to 
He i 1.083|um. 

The mass-loss rate is actually derived as the ratio 
M/a/7j where / is the volume filling factor that can be 
constrained by fits to the electron scattering wings of the 
helium line profiles (following Hillier 1991). We conclude 
/ ^0.1, and can definitely exclude homogeneous mass- 
loss in HD 50896. In addition, the mid-IR continuum 
slope also reacts to different filling- factors. 



TABLE 1 

Comparison of stellar parameters for EZ CMa (HD 50896, 

WN4b) derived here (labeled 'Model') with those 

determined previously by Schmutz (1997, S97), allowing 

for the clumped nature of the wind in each case with a 

volume filling factor of /. 
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3.2. Spectroscopic Results 

We compare our synthetic model with far-UV, near- 
UV, optical, near-IR and mid-IR spectrophotometry 
(dereddened by E(B-V)=0.10 mag) in Fig. Overall 
the agreement between the spectral features and contin- 
uum are excellent from 0.09-37/im, with few major ex- 
ceptions. For helium, all lines are reproduced better than 
20%, with the exception of a few mid-IR lines, namely 
Hen 3.091/mi (7-6), 13.12/mi (11-10) and 22.17 (jm (13- 
12). Since fine-structure lines are not included in our 
synthetic spectrum, it is apparent that [Ne in] is a non- 
negligible contributor to the 15.5/zm feature, the 25.9/im 
line is dominated by [O iv] and there is no obvious iden- 
tification for the lines at 20.7/mi or 36.6/Km (the latter is 
not [Nem] 36.01/im). 

The stellar parameters derived for the WN4b star are 
presented in Tabled We estimate elemental abundances 
of N/He=2xl0~ 3 and C/He=5xl0~ 5 by number. Hillier 
(1988) previously estimated N/He< 4 x 10~ 3 , N/C~14 
and 0/N<3 for EZ CMa, the only previous study to 
attempt CNO abundance determinations. 

The only recent study of EZ CMa allowing for blan- 
keting and wind clumping was by Schmutz (1997), whose 
results we also include in Tabled Schmutz (1997) ap- 
proached the incorporation of line blanketing in a dif- 
ferent manner from the present approach. A Monte 
Carlo sampling technique was adopted, allowing for the 
effect of a much more thorough line opacity at the ex- 
pense of full consistency in the radiative transfer prob- 
lem, via the use of approximate ionization and excita- 
tion for metal species. In addition, Schmutz adopted 
a He II Lya A303 photon loss mechanism with a par- 
ticular parameter. Such a mechanism is intrinsic to 
the present analysis without the requirement of a pa- 
rameterization, providing spectral lines adjacent to He n 
A303.78 are accounted for in the calculation. Conse- 
quently the following metal lines within ^100 kms -1 of 
Lya are included - Qui A303.70, 303.80, Fevi A303.70, 



303.80, Nivi A303.71, Mnvi A303.72, Cav A303.74, Crv 
A303.82, Crvi A303.84 and Fev A303.91. 

Agreement between the two approaches is reasonable, 
given the differences in approach and choice of diagnos- 
tics. In both solutions, stellar luminosities (and hence 
bolometric corrections) are significantly higher than pre- 
vious non-LTE studies (e.g. Hamann et al. 1995), al- 
though Schmutz (1997) obtained a rather higher bolo- 
metric luminosity, due to a combination of his photon 
loss mechanism and more complete line opacity, whilst 
clumping corrected mass-loss rates are much lower. Con- 
sequently, the derived momentum ratio is found to be 
Mvoo/(L/c) = 7, versus 30-70 in earlier studies. The 
agreement in mass-loss rate is consistent with the claimed 
precision of 0.1 dex achieved by Schmutz (1997), of rel- 
evance to the determination of elemental abundances to 
follow. 

3.3. Elemental Abundances: Neon and Oxygen 

We present the fine-structure lines of [Nem] 15.5/im 
and [OIV] A25.9/nn in Fig. including the predicted 
fractional contributions from [Ne in] and [O iv] (dotted 
lines) obtained from, respectively, the red line profile of 
[Neni] which does not overlap with Hen (15-13), and 
the non-LTE predicted He n (23-19) strength. The ob- 
served line fluxes of each blend is presented in Table |3 
Sources of atomic data are given in Dessart et al. (2000) 
for [Nem], and Hayes & Nussbaumer (1983) and Blum & 
Pradhan (1992) for [Oiv]. We have determined elemen- 
tal abundances for these ions using the numerical tech- 
niques introduced by Barlow et al. (1988), and adapted 
to account for a clumped wind by Dessart et al. (2000) 
and Morris et al. (2000). With regard to clumping, we 
admit that the volume filling factor in the outer wind 
may be considerably different than that derived from 
optical and near-IR recombination lines (Runacres & 
Owocki 2002). 

The numerical expression for the ion number fraction, 
7j, in a clumped medium (with volume filling factor /) 
is (in cgs units) 
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where D is the stellar distance, I u i is the line flux of the 
transition with energy hv u i between upper level u and 
lower level I, with transition probability A u i. j e (=1.008) 
and T e (=14000K) are the electron number density and 
temperature in the line-forming region, and the mean 
molecular weight is /i (=4.04) for EZ CMa, with ran the 
mass of the hydrogen atom. The integral part, F U (T), is 
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fu(N e ,T) 



d\og{N e ). 
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where f u is the fractional population of the upper level. 
Following this technique, we derive Ne 2+ /He = 
1.2xl0~ 4 and 3+ /He = 4.2xl0~ 5 by number, assuming 
a fractional contribution of 67% and 85% to the observed 
15.5/im and 25.9/im lines respectively. If we were to 
neglect the predicted Hen contribution to the observed 
spectral features, we would obtain Ne 2+ /He = 1.8xl0 -4 
and 3+ /He = 4.8xl0~ 5 from this method. Applying 
the same technique based on the Schmutz (1997) spec- 
troscopic results, we would derive slightly lower abun- 
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Fig. 2.— Comparison bwtween de-reddened far-UV (HUT), UV (IUE), optical (Torres-Dodgen & Massey 1988), ncar-IR (ISO) and mid-IR 
(Spitzer) spectrophotometry of EZ CMa with our synthetic model (dotted), corrected for interstellar atomic hydrogen with log(N(HI)=20.7 
cm~ 2 (Howarth & Phillips 1987). Note that the predictions from our non-LTE analysis are included here, such that specifically the 
fine-structure lines are not synthesiscd. 



dances of Ne 2 +/He = l.l-1.4xl0" 4 and 3+ /Hc = 3.5- 
4.1 xl(T 5 . 

Of course, one needs to account for potential contri- 
butions from unseen ionization stages. Our model atmo- 
sphere identifies the following dominant ions at the outer 



model radius of ^200i?*, where log(n e /cm 3 ) ~9, namely 
He+, C 3 +, N 3 +, 3 +, Ne 2 +, Mg 2 +, Si 4 +, P 4 +, S 4 +, Cl 4 +, 
Ar 4 - 5 +, Ca 4 +, Cr 4 +, Mn 5 +, Fe 4 ' 5 + and Ni 5 +. The dom- 
inant ionization state in WN stars at large radii is pre- 
dicted to be higher than for WC stars, despite similar 
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Fig. 3. — Mid-IR fine-structure lines comprising [Nc in] 15.55/im 
and [Oiv] 25. 9 fim observed in the Spitzer IRS observations of 
EZ CMa, together with the predicted fractional coontribution from 
the fine-structure lines (dotted) themselves (see text). 



O/He ~ 3+ /Hc = 4.2-4.8xlCT 5 . From our nitrogen 
abundance derived above, we find N/O=40-48 by num- 
ber. What is the expected oxygen abundance in a H-frec 
WN star? From the spectroscopic analysis above, we de- 
rive N/C=40 by number. This may be compared with 
recent predictions for (initially) rotating massive stars 
(Meynet & Maeder 2003). At the end of H-burning, the 
6O-12OM0 tracks for initial equatorial rotation velocities 
of 300 kms -1 predict N/C=40-45 by number, in perfect 
agreement. At this stage, N/O=30-50 by number, such 
that the measured value lies in precisely the predicted 
range. 



TABLE 2 
Observed mid-IR fine structure line intensities (units of 
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ERG CM 2 S X ) IN EZ CMa (uncorrected for 



potential blends with He II lines), and adopted atomic 

parameters, including statistical weights of the upper and 

lower levels, U) u AND UIi, transition probability, A ui , AND 

COLLISION strength U u i AT T e = 14000K. 
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[Neml 
[Oiv] 



3po 3po 
r l r 2 
2po _2po 
r l/2 r 3/2 



15.55 3 5 5.99xl(T 3 1.65 7.45xl0~ 13 
25.89 4 2 5.19xl0~ 4 2.52 1.49xl0" 13 



stellar temperatures, since the effects of metal coolants 
is far less (Hillier 1988, 1989). Are these predictions con- 
sistent for the physical conditions at even lower densities 
of log(n e /cm 3 ) ~5 where the fine-structure lines form? 

The absence of [S iv] 10.5/im is consistent with a dom- 
inant ionization stage of S 4+ for EZ CMa, whilst [Ne 11] 
12.8/im line is absent, so Nc + <C Nc 2+ , in agreement 
with expectation. The contribution of Nc 3+ is uncer- 
tain, although it has a rather high ionization potential 
(IP = 97eV), such that we assume Ne/He« Ne 2+ /Hc = 
1.2-1.8 xl0~ 4 , depending on the contribution by Hell. 
This is in reasonable agreement with the revised Solar 
Ne abundance (Asplund et al. 2004) of Ne/He=2xl0~ 4 , 
adjusted for the H-depleted atmosphere of an early WN 
star. 

For oxygen, we expect that 3+ is the dominant ion- 
ization stage, but are unable to verify the absence of fine- 
structure 2+ lines, which lie at 51.8 and 88.2/zm, long- 
ward of the IRS passband. Nevertheless, from the con- 
sistency with other relevant ions of Ne and S, we adopt 



4. SUMMARY 

We present Spitzer IRS spectroscopy for EZ CMa 
(HD 50896 WN4b), the first early WN star to be ob- 
served spectroscopically in the mid-IR. In addition to 
numerous stellar wind lines of He 11, fine structure lines 
of [Nelll] and [Oiv] are revealed, permitting an esti- 
mate of the asymptotic wind velocity, ~1860 kms -1 , 
plus elemental abundances. We carry our a spectro- 
scopic analysis of HD 50896, allowing for metal line 
blanketing and wind clumping, revealing generally ex- 
cellent agreement with the IRS spectroscopy, plus stel- 
lar parameters supporting the earlier study of Schmutz 
(1997), based on an alternative line blanketing tech- 
nique. We also derive N/He~ 2 x 10~ 3 by number, 
with N/C~40. From the derived clumped mass-loss rate, 
we obtain Ne 2+ /He=1.2xl0 -4 , allowing for the contri- 
bution by neighboring He n wind lines, in comparison 
with Ne/He=2xl0~ 4 for a severely H-depleted environ- 
ment, adjusted for the recent downward revision in the 
Solar neon abundance by Asplund et al. (2004). An 
indication of the uncertainty in the fine-structure de- 
rived Ne abundance may be obtained by adopting the 
spectroscopic results for EZ CMa from Schmutz (1997), 
leading to a 10-20% difference. In addition, we obtain 
O 3+ /He=4.2xl0~ 5 , i.e. N/0^48, in excellent agreement 
with evolutionary predictions for massive stars at the end 
of core H-burning. 
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